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The physical limit with which a cell senses external ligand concentration corresponds to the perfect 
absorber, where all ligand particles are absorbed and overcounting of same ligand particles does 
not occur. Here we analyze how the lateral diffusion of receptors on the cell membrane affects the 
accuracy of sensing ligand concentration. Specifically, we connect our modeling to neurotransmission 
in neural synapses where the diffusion of glutamate receptors is already known to refresh synaptic 
connections. We find that receptor diffusion indeed increases the accuracy of sensing for both the 
glutamate AMPA and NDMA receptors, although the NMDA receptor is overall much noiser. We 
propose that the difference in accuracy of sensing of the two receptors can be linked to their different 
roles in neurotransmission. Specifically, the high accuracy in sensing glutamate is essential for the 
AMPA receptor to start membrane depolarization, while the NMDA receptor is believed to work in 
a second stage as a coincidence detector, involved in long-term potentiation and memory. 
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INTRODUCTION 



Biological cells live in very noisy environments from 
which they receive many different stimuli. The survival 
of a cell highly depends on its ability to respond to such 
stimuli and to adapt to changes in the environment. A 
fundamental role in this task is played by membrane re- 
ceptors, proteins on the cell surface which are able to 
bind external ligand molecules and trigger signal trans- 
duction pathways. Although the precision with which 
a cell can measure the concentration of a specific lig- 
and is negatively affected by many sources of noise , 
several examples exist in which such measurements are 
performed with surprisingly high accuracy. In bacterial 
chemotaxis for instance, the bacterium Escherichia coli 
can respond to changes in concentration as low as 3.2 nM 
Q, corresponding to only three molecules in the volume 
of a cell. High accuracy is observed also in spatial sensing 
by single cell eukaryotic organisms as e.g. in the slime 
mold Dictyostelium discoideum, which is able to sense a 
concentration difference of only 1 — 5% across the cell 
diameter Q, and in Saccharomyces cerevisiae (budding 
yeast), which is able to orient growth in a gradient of a- 
pheromone mating factor down to estimated 1% receptor 
occupancy difference across the cell [7|. Spatial sensing 
is also efficiently performed by growing neurons, lympho- 
cytes, neutrophils and other cells of the immune system. 




FIG. 1: (Color online) Synaptic junction during neurotrans- 
mission. A vesicle containing a neurotransmitter (e.g. glu- 
tamate for excitatory synapses), releases its content into the 
synaptic cleft. Receptors on the post-synaptic surface bind 
the neurotransmitter, increasing the likelihood of the propa- 
gation of the electric impulse (action potential). Red arrow 
indicates the incoming potential, yellow arrow the set of reac- 
tions leading to propagation of potential, small green arrows 
indicate receptor mobility. 



There has been substantial progress in our theoretical 
understanding of the accuracy of concentration sensing. 
In 1977, Berg and Purcell were the first to point out 
that ligand sensing is limited by ligand diffusion, pro- 
ducing a "counting noise" at the receptors Bialek 
and Setayeshgar Q (and later others [10|]) used the Fluc- 
tuation Dissipation Theorem (FDT) to separate the con- 
tributions from random binding/unbinding events from 
rebinding due to diffusion of ligand molecules. Subse- 
quently, Wingreen and Endres showed that the perfect 



absorber is the absolute fundamental limit of ligand sens- 
ing accuracy; it does not suffer from ligand unbinding and 
rebinding noise [l(| EH • 

In previous work we analyzed the contribution of endo- 
cytosis, i.e. the internalization of cell-surface receptors, 
to the accuracy of external ligand concentration sensing 
[IH ■ We showed with a simple model that internalization, 
by making the cell act as an absorber of ligand, increases 
sensing accuracy by reducing the noise from rebinding 
of already measured ligand molecules. In this paper we 
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consider the effect of iateral diffusion of receptors on the 
accuracy of ligand concentration sensing. 

In a situation in which ligand sensing is concentrated 
in a specific region on the cell membrane (signaling 
hotspot), receptors bind ligand molecules locally but may 
release them remotely from the region of interest, thus 
preventing those particles from rebinding locally. Noise 
from overcounting could therefore be reduced, potentially 
increasing the accuracy of sensing. 

A potential implementation of such localized sensing 
occurs in neural synapses, responsible for transmission 
of action potentials and short-term synaptic plasticity 
(see Fig. 1 for a schematic description). At such crit- 
ical places in the central nervous system, the accuracy 
of sensing ought to be important [l3| . At a synapse the 
propagation of an action potential arriving from the pre- 
synaptic neuron causes the opening of ion channels on 
the pre-synaptic membrane which increases calcium-ion 



inflow 14]. Calcium ions trigger a biochemical cascade 



which results in the formation of vesicles containing neu- 
rotransmitter, e.g. glutamate, in the case of excitatory 
synapses. One or more vesicles eventually fuse with the 
membrane at a point on the pre-synaptic surface to form 
a fusion pore. Glutamate is subsequently released and 
diffuses into the synaptic cleft, reaching the receptors on 
the post-synaptic surface on the opposite side of the cleft. 
These receptors, mainly the AMPA and NMDA receptors 
(named after the agonists a-amino-3-hydroxy-5-methyl- 
4-isoxazolepropionic acid and N-methyl-D-aspartic acid 
respectively), cause ion channels on the post-synaptic 
membrane to open, facilitating the propagation of the 
action potential. 

The AMPA receptor (AMPAR) has a tetrameric struc- 
ture with four binding sites for glutamate. Binding glu- 
tamate on two sites leads to conformational change and 
opening of a pore. More occupied binding sites imply a 
higher current through the channel, mainly of calcium, 
sodium and potassium ions. The NMDA receptor (NM- 
DAR) has a tetrameric structure as well, but its chan- 
nel needs coincidental binding of glutamate and glycine 
molecules on dedicated sites for it to open. A small mem- 
brane depolarization is furthermore necessary to clear out 
the magnesium ions blocking the channel. This third 
condition makes the NMDA receptor a type of coinci- 
dence detector for membrane depolarization and synap- 
tic transmission, playing an important role in memory 
and learning. Interestingly, AMPA receptors diffuse on 
the post-synaptic surface unusually fast, especially when 
ligand is bound 15]. In addition to refreshing the synap- 
tic plasticity [16(, such diffusion may also increase the 
accuracy of sensing in neurotransmission. 

This paper is organized as follows: in Sec. II we re- 
view the case of a single immobile receptor; in Sec. Ill we 
consider receptor diffusion on the membrane and derive 
the dynamical equations determining the spatial concen- 
tration of ligand and the occupancy of the receptor. In 
Sec. IV we derive the stationary solution for the receptor 
occupancy, and in Sec. V, using a non-equilibrium ap- 



proach based on the effective temperature, we derive the 
accuracy of sensing. In Section VI we describe how our 
theoretical results connect to neurotransmission in neu- 
ral synapses. The final section is devoted to an overall 
discussion including cases of further biological relevance. 
Technical details are provided in two appendices. 



II. REVIEW OF THE SINGLE RECEPTOR 

In this section we review previous results for a single, 
immobile receptor. As depicted in Fig. 2, such a recep- 
tor can bind and release ligand with rates k + c and k-, 
respectively. The kinetics for the occupancy n(t) of the 
receptor are therefore given by 



dn(t) 
dt 



= k + c[l - n(t)] - k_n(t), 



(1) 



where the concentration of ligand, c, is assumed uniform 
and constant. The steady-state solution for the receptor 
occupancy is given by 



(2) 



K 



D 



with Kr> = k-/k+ the ligand dissociation constant. The 
rates of binding and unbinding are related to the (nega- 
tive) free energy F of binding through the detailed bal- 
ance 



(3) 



with T the temperature in energy units. In the limit of 
very fast ligand diffusion, i.e. when a ligand molecule is 
immediately removed from the receptor after unbinding, 
the receptor dynamics are effectively decoupled from the 
diffusion of ligand molecules and hence, diffusion does 
not need to be included explicitly. 

Following Bialek and Setayeshgar [§], the accuracy of 
sensing is obtained by applying the FDT [l7[, which re- 
lates the spectrum of fluctuations in occupancy to the 
linear response of the receptor occupancy to a pertur- 
bation in the receptor binding energy. Furthermore, at 
equilibrium the fluctuations in occupancy can be directly 
related to the uncertainty in ligand concentration using 
Eq. @. After time-averaging over a duration r much 
larger than the correlation time of the binding and un- 
binding eve nts, the normalized uncertainty of sensing is 
given by H [l| 



1 



A;_|_c(l — n)T 2ttD^,cst ' 



(4) 



where the right-hand side is obtained for diffusion-limited 
binding [lj|], i.e. when k + c(l — n) — > AncD^s, with D% 
the diffusion constant and s the dimension of the (spher- 
ical) receptor. Eq. (0| shows that the uncertainty is 
limited by the random binding and unbinding of ligand. 
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FIG. 2: (Color online) Single receptor, immobile at position 
xo, binds and unbinds ligand with rates k+c and respec- 
tively. 

The accuracy of sensing is defined as the inverse of the 
uncertainty. 

In the case where diffusion of ligand is slow, ligand 
binding to the receptor is perturbed [9]. The kinetics 
of receptor occupancy and ligand concentration are de- 
scribed by 

dn(t) 



dt 
dc(x, t) 



k+c(xo, t)[l — n(t)] 



D 3 V 2 c(x,t) - 5(x-x ) 



fc_n(f) 
dn(t) 



(5a) 
(5b) 



dt K ' ' ' dt 

where xq indicates the position of the receptor. In the 
last term in the second equation, the Dirac delta func- 
tion S(x — xq) describes a sink or source of ligand at x 
corresponding to ligand binding or unbinding from the 
receptor, respectively. 

Following a similar procedure as in the previous case, 
the uncertainty of sensing is given by nm 

((Sc) 2 )r 2 



fc+c(l — n)r 
3 



1 

itD^cst 



(6a) 
(6b) 



2ttD 3 cst 

where the first term on the right-hand side of Eq. (paj) 
is the same as in Eq. (j4j, while the second term is the 
increase in uncertainty due to diffusion. This term ac- 
counts for the additional measurement uncertainty from 
rebinding of previously bound ligand to the rece pto r. For 
diffusion- limited binding, one obtains Eq. (fob)) 

Comparison of Eqs. (jU and © shows that removal of 
previously bound ligand by fast diffusion increases the 
accuracy of sensing, since the same ligand molecule is 
never measured more than once. 

The uncertainty in Eq. Q can be further reduced 
by a factor of two, consequently increasing the accuracy 
of sensing by the same factor, by considering the fun- 
damental limit. In this limit, ligand particles are ab- 
sorbed, and hence only binding noise contributes. For 
diffusion-limited binding, the fundamental limit is given 

by idaim 



{(5c) 



I 




FIG. 3: (Color online) Ligand-receptor binding by mobile re- 
ceptors characterized by diffusion coefficient D2 (green ar- 
rows). Ligand (red dots) is concentrated only in a small area 
AA around xo due to e.g. local release of giutamate from vesi- 
cle fusion pore with rate koo . Ligand can bind to receptors at 
xo, but unbind everywhere on the membrane. Parameter D3 
describes the diffusion coefficient of the ligand. 



calculated from the diffusive flux to a small sphere of 
radius s, which represents the receptor. 



III. RECEPTOR DIFFUSION MODEL 

In this model, depicted in Fig. 3, receptors are free 
to diffuse on the 2-dimensional (2D) membrane surface 
with diffusion coefficient D^. A receptor located in the 
limited area A A centered at xq, can bind and unbind 
ligand of average concentration c with rates fc+c and 
respectively. When a receptor leaves this area and dif- 
fuses to the remaining large region of negligible ligand 
concentration, it can either release ligand or diffuse back 
in the area. 

We consider a situation in which receptor diffusion has 
reached a stationary condition resulting in a uniform den- 
sity po of receptors on the membrane, and this condition 
is not changed by interaction with ligand molecules. Lig- 
and binding and unbinding only affects the number of 
occupied and unoccupied receptors with respective den- 
sities pb and p u , given by 



Po = Pb(r) + Pu(r), 



(8) 



4ttDscst 



(7) 



where r = (2, y) is a position on the 2D membrane sur- 
face. For purposes of comparison with the single immo- 
bile receptor model, we also assume that there is always 
exactly one receptor in the area A A around xq — (tq, Zq), 
i.e. po = l/A. 

The dynamics of such a system is captured by the fol- 
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lowing set of coupled differential equations 
dp b (r,t) 



dt 



D 2 V 2 p b {f, t) - k^pbif, t)+ (9a) 
k + [p - p b (r, t)} c(r, z , t) AA <5(r - fo) 



dc(x, t) 
dt 



D 3 V 2 c(x,t) - 5(z - z Q ) 



dp b (r,t) 
dt 



-D 2 V 2 p b (f,t) 



(9b) 



where rate constant fcoo ensures the replenishment of lig- 
and molecules, x — (x,y,z) — (f,z) indicates a location 
in 3D space. Specifically, p b (f,t) is the density of bound 
receptors and c(x, t) is the ligand concentration, both de- 
pending on space and time. An additional equation for 
p u (f,t) is not necessary due to receptor conservation in 
Eq. (8). 

After normalization by the average density, using 
n(f,t) = p b (f, t)/po in Eq. ©, we obtain 



dn(f, t) 

dt 
dc(x, t) 

dt 



D 2 V 2 n(f,t) +k+[l- n(f,t)} c(f, z = z ,t) AAS(r- f Q ) - k-n{r,t) 



= D 3 V 2 c(f , t) - p S(z - z ) 



dn(f, t) 
dt 



-D 2 V 2 n(r,t) 



(10a) 
(10b) 



These equations describe the coupled dynamics of 
the diffusing receptors and their occupancy with ligand. 
In summary, receptors can bind ligand molecules only 
within area AA but can release them anywhere on the 
membrane. The solution is provided in the following sec- 
tions. 



IV. NON-EQUILIBRIUM STATIONARY 
SOLUTION FOR RECEPTOR OCCUPANCY 

It is possible to extract an exact analytical expression 
for the stationary solution n(r) for the receptor occu- 
pancy from Eq. (|L0a|) . Setting the right-hand side to 
zero and Fourier Transforming both sides leads to 



n(q) = 



k + [l-n(r )]c(x a )AA 
D 2 q 2 + k- 



(11) 



where q is the wave vector of magnitude q = \q\. The 
inverse Fourier Transform, upon introducing a cut-off 
A ~ l/s due to the receptor size, leads to the follow- 
ing solution 



n(r) = k + [l - n(fo)]c(x ) 



AA 

4ttD ? 



log 1 



D 2 A< 



(12) 



This stationary solution corresponds to a non- 
equilibrium condition since ligand is continuously taken 
away from area A A. Setting A A = Att/A 2 ~ s 2 and 
evaluating Eq. (fT2"|) in fo leads therefore to the following 
value for the occupancy at fo 



1 



k+c(x 0/ 



AA 



log (l 



4ttD 2 
AAk- 



(13) 



which is plotted in Fig. 4 and, in the limit D 2 —> 0, leads 
to the recovery of the immobile receptor solution Eq. ^ , 
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FIG. 4: (Color online) Receptor occupancy at fo- Eq. fj 13p 
for n is plotted as a function of lateral diffusion coefficient D2 
and unbinding rate fc_ for a range of experimentally relevant 
values in neurotransmission [top surface). Specifically, the 
k- values for AMPA and NMDA receptors are indicated (see 
Table I). Corresponding plot for n(l — n) is also shown (bottom 
surface). 



namely uq — k + c/(k^ + k + c). In order to quantify the 
deviation from equilibrium we apply Fick's first law to 
calculate the flux J of occupied receptors 



J = D 2 Vn(f,t) 



(14) 



out of region AA [19| . Integrating Eq. (|14l) over the bor- 
der of AA allows us to define an effective "internalization 
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rate" for removal of ligand in analogy to [12J , via 
D 2 



kfn 



AA 



\7n(f)dr 



' aa 



V 2 n{f)d 2 r 



AA 



= — / d 2 f[k+c(l - n)AAS(r) - k-fi] , (15) 
J AA 

where the second equality follows from integration by 
parts and the last equality results from the stationary 
condition for Eq. (|10ap . Carrying out the last integration 
leads to the following final expression for the effective 
internalization rate 



fc + c 



l-n 



47rD, 



AA log 1 



(16) 



where we used Eq. (fl3|) to express n in terms of the re- 
ceptor diffusion coefficient. Using condition A A — 4ir/A 2 
we can define a modified unbinding rate 



fe+c(l — n) 



log (l + ^ 



(17) 



and extend the equilibrium condition Eq. (3) to the sta- 
tionary non-equilibrium condition via an effective tem- 



J 



perature T e defined by [U 



k+c f 
— — = e T = . 



(18) 



This effective temperature allows us to generalize the or- 
dinary FDT to a stationary non-equilibrium condition 
and to derive the receptor accuracy of sensing in the next 
section. 



V. ACCURACY OF SENSING 

Here we we derive the accuracy of sensing for diffus- 
ing receptors. Starting from Eq. (|10p we consider an 
expansion to first order around the stationary solutions 



n(f, t) = n + Sn(f, t) 
c(x, t) — c + 5c(x, t), 



(19a) 
(19b) 



where n and c are the stationary solutions of Eq. (|10l) , 
which, in this case, have a spatial dependence, i.e. n = 
n(f) and c = c(x). Hence, the linearized equation for the 
receptor occupancy and ligand concentration are given 
by 



d[Sn(r,t)} 
dt 

d[Sc(x,t)} 
dt 



D 2 V 2 6n{r,t) 



_ 6k+(t) _5k-(t) 
(1 — n)oc(r, zq, t) — con(r, t) + (1 — n)c — : n — : 



D 3 \7 2 Sc(x,t) - 6(z- z ) Po < 



k+AA8{r- fh) - k-6n(r,t) 
(20a) 

(20b) 



r 



where we assumed that fluctuations in the bind- binding free energy, resulting in 
ing/unbinding rate constants occur as well. This math- 
ematical trick allows us to introduce fluctuations in the 



receptor-binding free energy and so to apply the FDT d[Sn(f t)] 
0, E3] • Using the stationary non-equilibrium condition 
Eq. (fT5)l . we obtain the fluctuations in the rates and 
binding free energy 



dt 



D 2 V 2 6n(r,t) - k-5n(r,t) + 



SF 

fc+c(l - n) — 



+fc + (l — n) Sc(f, Zq, t) — k+c5n{f, t) 



AAS(r- r"o ). 



(22) 



6k+ 5k_ _ SF 

fe-L K_ Tp 



(21) 



Fourier Transforming Eq. (|22|) and setting q = we 
solve for 8h(q,u>) and obtain 



Using this equation allows us to replace fluctuations in 
the rate constants in Eq. (I20a[) with fluctuations in the 



5h(q,u) = G(uj,r ,z ) 



D 2 q 2 



k— — iu 



(23) 
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where, for convenience of calculation, we have defined the 
following function 



G(u,r ,z ) 



fc + (l — n) Sc(f , z ,lu) — k + cSn(f ,ui)+ 
6F(lo) 



(1 - n)c 



T, 



AA. 



(24) 



Inverse Fourier Transforming Eq. (|23l) leads to the fol- 
lowing expression for the spectrum of the fluctuations in 
occupancy 



5h(f ,oj) 



dq 



(2^ 
= G(uj,r ,z Q ) 



lq - ro 6n(q,uj) 
dq 



(25) 



1 



(2tt) 2 D 2 q 2 +k-- iuj 
G(w,ro,zo)Si(w), 



where we defined 
Ei(w) 



dq 



1 



(2tt) 2 D 2 q 2 + k- - iu) 



(26) 



with the integral provided in Appendix A. In order to 
remove the ligand concentration in Eq. (1241) . we Fourier 
Transform Eq. (|20bjl . lead ing to 

6c(q,uj) = po n ™~ D 2 2 f ■ 6n(q,Lj)e^ z °. (27) 

Inserting Eq. ([2^]) in Eq. (|27p and inverse Fourier Trans- 
forming leads to the following expression for the spectrum 
of the fluctuations in ligand concentration 



6c(f ,z ,uj) = I 



-iq±z 



dq 

2tt ~ J (2tt) 2 

G(w,r , z )p T, 2 (ui) 



e - i9 "°(5c(g,w) 
(28) 



with 



S zH = / — 



2^ J (2tt) 2 £ 3 (<Z 2 + <£) - iw g 2 + *p - ^ 

2 (29) 

provided in Appendix B. Inserting Eqs. (|25|) and (|28|) 
in Eq. (|2"41 leads to a closed equation for the function 

G(uj,r Q ,z ) 



G(uj,r ,z ) 

— fc_cEi (w)G(w, fo, ^o) + fe+c(l — n) 



fc+(l - n) ( o E 2 (w)G(a;,ro,zo)-l- (30) 



which can be solved for G(ui, fo, 2rj)- This leads to 
k+c(l-n)AA 8F(oj)/T e 



G(u,r ,z ) 



1 - fc+(l - n)S 2 (w) + fc + cEi(w)AA : 



(31) 



g .3 2 - 
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FIG. 5: Receptor-occupancy fluctuations ((Sn) 2 ) T in units of 
t as a function of rj = 16D 2 /(s 2 fc_) for AMPAR (a) and 
NMDAR (b). Shown are (solid line) plot of {{Sn) 2 ) T as given 
by Eqs. (JMJ and po). (dashed line) plot of Eq. using the 
effective internalization rate Eq. (|16|) . and (vertical thin line) 
experimental value of r\ for each receptor type (see Table I 
for all parameters value) . Glutamate concentration was set to 
c=0.1mM, here and in Fig. 6. In order to produce a significant 
effect from ligand rebinding, we used D3 = Q.lp,m 2 / s here and 
in subsequent plots, as may result from binding of glutamate 
to neuroligins, neurexins and other cleft proteins. 



where condition po = 1/AA was applied. Using Eqs. ([23] 
and (|3ip . a final expression for the response is obtained 



(32) 



5F(u) 



6F(u) 



k+c(l - n)AiEi(w) 



T e 1 - fc+(l - n)£ 2 M + fc+cSi(cj)AA 



Applying the generalized FDT we derive the spectrum 
of the fluctuations in receptor occupancy 5n(ro,uj) from 
the response, namely 



2kT P 



Tm 



Sn(r ,u>) 
SF(u>) 



(33) 



By replacing the expressions for Ei(oj) and ~S 2 (uj) in Eq. 
), taking the imaginary part and setting w ~ 0, we 
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obtain for the zero- frequency limit of the power spectrum 
fc+c(l - n)AAA 2 



5,(0) 



2irk 2 _ (1 + 77) 



1 



fc + (l-n)A 



1 fc + (l-n)A / U 3 A 2 - r?fc- ^ 



16irD 3 



(1 - 4,) 
1 - A 



(34) 



where we conveniently defined the following functions of 
the dimensionless parameter r\ = D2A 2 



log (1 + 77) 



A... 



n 



arctan Jr\ 



(35) 



Note that in Eq. (|34j) . occupancy n also depends on i] 
via Eq. (fl3|) . The zero-frequency spectrum is related 
to the receptor occupancy fluctuations, averaged over a 
time t ^> kZ , (/c+c) -1 , by the following relation: 



((5n)\ ~ 5,(0)/ 7 



(36) 



In Fig. 5 we plot the occupancy fluctuations, given by 
Eq. p6|) as a function of i), i.e. the effective diffusion 
strength. This is done for the two main receptors in- 
volved in neurotransmission, namely AMPAR and NM- 
DAR, with parameters taken from Table I. For AMPAR, 
the effect of diffusion is that of decreasing the fluctuations 
in occupancy due to reduced overcounting of previously 
bound molecules (Fig. 5a). For small unbinding rates 
k- though, such as occurs for the NMDA receptor, an 
actual increase in the noise is observed for a range of 
physiologically relevant diffusion strengths. 



In order to make the origin of this effect more transpar- 
ent, we also plot the occupancy fluctuations in Fi g. 5 as 
calculated for the case of receptor internalization [l2| 



((Sn) 



2n 2 (l-rl) 



(37) 



where n — k + c/(k + c + fc_ + kf), with k\ defined in Eq. 
(TIT))) . The experimental unbinding rate of NMDAR leads 
to values for n such that starting from D2 = 0, the func- 
tion ri(l — n) increases, reaches a maximum and then 
decreases again. Such a maximum is avoided by AM- 
PAR whose spectrum monotonically decreases when D2 
increases (see Fig. 4). The comparison between recep- 
tor diffusion and receptor internalization in Fig. 5 shows 
qualitatively similar results, confirming the effective lig- 
and removal by receptor diffusion. 



A. Equilibrium limit for vanishing D2 

Here we show that Eq. (|34l) leads to the expression 
for the single receptor obtained in Refs. 0, [l8[ in the 
limit for D2 — > 0, i.e. r\ — > 0. Taking the limit and using 
A A = 4ir/A 2 , we indeed obtain 



hm^(0)= (fc _ + fc+5)2 



fc + (l-n )A 



(38) 



which is identical to the result for the single receptor for 
A = 4/s. 



Da 
k+ 
k- 



GLUTAMATE 



300 fan /s(acqueous solution) [20] 
10/im 2 /s (synaptic cleft) [21, 22] 
0.1-1 mM [20, 23, 24J 
1 ms [25] 
8 nm [14, 26] 



AMPAR 



0.028-0. l/jm'/s [27, 28] 
4-10 xlO^-V 1 ^, 29] 



2-8 x 10'' 



[20, 29] 



NMDAR 



0.021/im7s [27, 28} 
5-8.4 xlO^-^ -1 ^, 30] 



[29, 301 



TABLE I: Summary of experimentally determined parameters 
for glutamate ligand (top), AMPA receptor (bottom left) and 
NMDA receptor (bottom right). 



B. Near-equilibrium result for small D2 

We derive here the exact result, Eq. (|3"4"|) . in the limit 
of slow receptor diffusion on the membrane, using r\ = 
D2/{s 2 k-) as a small parameter. Specifically, we are 
interested in the accuracy of sensing, which is derived 
from the normalized variance for the ligand concentration 
of ligand as in [§] 



((Sc)\ ((Sn) 2 ) T 



1 



S v (0) 



n 2 (l — n) 2 n 2 (l — n) 2 



(39) 



where the right-hand term follows from Eq. (|36j) . Ex- 
pansion to first order in parameter 77 leads to 



k + c(l — 71q)t 



(1 - 8n ?7) + 



1 



1 - 

r 



16k 2 n 2 
„ 7 - 2n ° +8"o- 



3fc_ — k+c fc_ 



H s V 

3fc + c 2D 3 1 ' 



(40) 



Comparison with Eq. (6a) shows that the correction der: (i) the first term is reduced due to a reduction in n, 
due to lateral receptor diffusion is twofold to first or- 
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FIG. 6: Uncertainty ((Sc) 2 ) T /c 2 as a function of r\ = 
16D 2 /(s 2 fc_) for the AMPA receptor (a) and the NMDA re- 
ceptor (b) in units of r. Shown are (solid line) plot of Eq. 
f!3Q P as obtained from the exact solution Eq. (I34|l . (dashed 
line) linear approximation for small r\, and (dotted line) con- 
tribution to the accuracy of sensing from only the first term in 
Eq. (|40p . Also shown are (vertical thin line) the experimen- 
tal value of r\ for each receptor (see Table I) and (horizontal 
dot-dashed line) the fundamental physical limit Eq. 



(ii) the second term is reduced for most parameter val- 
ues due to a reduction in rebinding of already measured 
ligand molecules. 

In Fig. 6 we plot the uncertainty in sensing for both 
AMPAR and NMDAR. In both cases the total effect of 
lateral diffusion is that of decreasing the uncertainty in 
the concentration and therefore of increasing the accu- 
racy of sensing. Receptor diffusion for AMPAR (Fig. 
6a) mainly reduces the rebinding term, while for NM- 
DAR (Fig. 6b) the first term from binding and unbind- 
ing is reduced. In the latter case, the maximum in the 
occupancy fluctuations (see Fig. 5) is removed through 
fi 2 (i — n) 2 in the denominator of Eq. For large val- 

ues of r), the occupancy fluctuations and hence the accu- 
racy of sensing become unphysical (zero, i.e. below the 
physical limit, or diverge) as the effective temperature 
breaks down far from equilibrium. (For large rj, recep- 
tors remove ligand molecules efficiently, which introduces 
large non-equilibrium ligand fluxes. These cannot be rep- 




20 30 
current I (pA) 

FIG. 7: Quantal transmission, (a) Histogram of experimen- 
tal EPSC (excitatory post-synaptic current) data taken from 
Ref. [Hj] and fitted by three Gaussian distributions, (b) Dis- 
tribution of predicted EPSCs. Shown are individual Gaus- 
sian distributions corresponding to one, two and three vesi- 
cles (from left to right; thin lines), as well as envelope func- 
tions (sums with equal weighting; thick lines). Calculation is 
done for AMPA receptors using Eqs. (HU gl} , Eq. d39j for 
the variances of diffusing (thin solid lines) and Eq. (6a) for 
immobile (thin dashed lines) receptors, assuming giutamate 
concentration c n = ct + nc v for n = 1,2,3 vesicles, back- 
ground cj, = 0.02mM, vesicle concentration c v — 0.04mA/ 
and N = 100 receptors [ii liEll . Furthermore, Eq. (fiTj) with 
parameters from Refs. 34, 35] was used. 



resented by an effective temperature.) 

In summary, lateral receptor diffusion affects the accu- 
racy of sensing in two ways: (i) by reducing the station- 
ary value for the receptor occupancy, and (ii) by reduc- 
ing the noise from rebinding of already measured ligand 
molecules. 



VI. QUANTAL TRANSMISSION 



A remarkable property of synaptic transmission is that 
the amplitude of the synaptic response, i. e. the excita- 
tory post-synaptic current (EPSC), varies by an integral 
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multiple of a quantum. This property, first unveiled by 
the pioneering work of Katz [31| , was traced to the forma- 
tion of synaptic vesicles and subsequent pore formation. 
The latter leads to a release of glutamate into the synap- 
tic cleft [321 ] . This quantization is considered a way to 
minimize the variability of the response from each con- 
tributing synapse and therefore to allow efficient neural 
computation. This variability is quantified by the coeffi- 
cient of variation (CV) defined as the standard deviation 
devided by the mean value. The reason for a large CV 
is controversial (cf. Fig. 7a), but can be attributed to a 
variation in the size of the vesicles carrying glutamate, or 
to the fact that large EPSCs are generated by multiple 
vesicles [IH Hf|. The experimental detection of EPSCs 
with an amplitude distribution of more than one peak is 
consistent with the latter explanation [33| . 

In Fig. 7a, the experimental data taken from Ref. [34j 
show EPSC peaks as a result of repetitive suprathresh- 
old stimulation of granule cell somata in CA3 pyramidal 
cells. Some stimuli did not evoke an EPSC (termed fail- 
ures). The histogram of successful events can be fitted 
with three Gaussian distributions, showing the quantal 
nature of the EPSC. The peaks of the Gaussians corre- 
spond to the discrete multiple values 7, 14, and 21 pA, 
respectively. 

We can compare this distribution to the result pre- 
dicted by the accuracy of sensing. Specifically, we would 
like to know if fast diffusion helps resolve the EPSC 
peaks. A simple model based on the structural properties 
of AMPAR was previously used to derive the dynamics 
of glutamate binding and unbinding, as well as channel 
opening and closing (34|. Based on this model the fol- 
lowing relation connecting the current through the ion 
channel of AMPAR and the concentration of glutamate 
was deduced: 



I = Io 



1 



1 + (A/c)' 



(41) 



where Iq, A, and the Hill coefficient n are parameters 
provided in Ref. [34|,|35j]- Based on Eq. (|4"Tj) the following 
relation connecting the width in the distribution of EPSC 
currents to the fluctuations in glutamate concentration 
can be derived: 



SI = IqhX 



(x/ c y 



[i + {xi c y 



tSc. 



(42) 



We further assume that a number N of receptors inde- 
pendently contribute to the measurement, i. e. 



((Sc) 2 ) T . N = ((Sc) 2 ) T 
c 2 Nc 2 ' 



(43) 



with the uncertainty of a single receptor given by Eq. 
(f3"9"| . Averaging over N, therefore, further reduces the 
uncertainty in Eq. (|43[) and hence increases the accuracy 
of sensing. 

Figure 7b shows that the distributions of predicted EP- 
SCs is in qualitative agreement with the experiments. 



From our calculations we found that the effect of lateral 
receptor diffusion is to increase the resolution of peaks 
and so also the ability of the synapse to count the num- 
ber of vesicles released into the synaptic cleft. Note that 
in order to resolve the EPSC peaks experimentally, the 
probability of vesicle release from the presynaptic side 
was drastically reduced. This was achieved by modulat- 
ing the ratio of Ca + /Mg + ions in the extracellular so- 
lution, resulting in a 50-60% reduction of peak currents 



VII. DISCUSSION AND CONCLUSIONS 

Our previous work on the fundamenal physical limit 
of ligand sensing led us to conclude that receptor 
diffusion, by removing bound ligand from a region of 
interest, may reduce the local overcounting of ligand 
molecules, therefore potentially increasing the local ac- 
curacy of sensing. In this paper we set out to investigate 
this possibility by constructing a mathematical model, 
that contains all the necessary ingredients to describe 
the role of receptor diffusion in ligand sensing. We have 
considered a receptor, that can bind and unbind ligand 
molecules as well as diffuse on a 2D membrane. Lig- 
and is allowed to diffuse in 3D space. Using this model 
we derived the fluctuations in receptor occupancy in the 
area of interest via a generalization of the FDT with the 
introduction of an effective temperature [Hj]. This ulti- 
mately allowed us to derive an equation for the accuracy 
of sensing. 

We applied our model to the biologically relevant case 
of glutamate receptors, which are responsible for trans- 
mitting action potentials at neural synapses. We found 
that for AMPAR and NMDAR, our model shows that the 
lateral diffusion of receptors increases their accuracy of 
sensing and hence may allow synapses to count the num- 
ber of released vesicles. However, the local occupancy 
fluctuations are decreased only for AMPAR; for NMDAR 
these fluctuations are increased. This difference is due to 
the smaller unbinding rate of the latter receptor. Conse- 
quently, the accuracy of sensing for NMDAR is an order 
of magnitude smaller than for AMPAR. It is important 
to remark that these results are consistent with the differ- 
ent roles these receptors play. In fact, the transmission 
of an excitatory potential through a synapse occurs in 
several steps. Initially, it is AMPAR which senses gluta- 
mate and, by letting in ions, starts a small depolarization 
of the membrane. This in turn allows NMDAR to open 
its channel and start an even bigger influx of ions (es- 
pecially Ca + ). This stimulates the production of more 
AMPARs and so increases the strength and sensitivity 
of the synapses. Sensitivity to glutamate is important 
for AMPAR, while the roles of NMDAR are coincidence 
detection and amplification, important for long-term po- 
tentiation and memory. 

In order to derive the accuracy of sensing for diffusing 
receptors, we made a number of symplifying assumptions. 
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We neglected fluctuations in the total receptor density 
and vesicle size, assumed to be constant [TH H^j. Fur- 
thermore, we assumed a constant ligand concentration 
in the hotspot area AA during the receptor measure- 
ment time, while in reality the concentration profile will 
broaden out. Additionally, we made the simplifying as- 
sumption that multiple vesicles are released in the same 
spot, so that released glutamate concentrations add up, 
and that other effects such as spillover from neighboring 
cells can be neglected [H, [43j]. Finally, we introduced 
the effective temperature Tl to gen eralize the FDT to 
non-equilibrum processes [12L l36l439l | . The approximate 
nature of T e , for which we neglected any potential time 
or frequency dependence, limits the quantitative valid- 
ity of our model to small deviations from equilibrium. 
In fact, for large receptor diffusion coefficients, the occu- 
pancy fluctuations and hence the accuracy of sensing be- 
come unphysical. Nevertheless, the mapping of receptor 
diffusion to an effective internalization process, already 
studied in Ref. [ljj], provides confidence in our method. 

In conclusion, in this paper we highlighted the role 
of diffusion in increasing the accuracy of sensing. This 
might be important for synaptic counting of glutamate 
vesicles. This hypothesis can be tested experimentally by 
reducing the mobility of receptors. This can be achieved 
by cross-linking or addition of cholesterol. Alternatively, 
the glutamate diffusion costant can be reduced by adding 
dextran. Such perturbations should lead to a broadening 
of the EPSC peaks, although signaling may be affected 
as well. Similarly to AMPAR [15], an increase in re- 
ceptor diffusion upon ligand binding was recently also 
observed in the pse udopods of migrating Dictyostelium 
discoideum cells [40j. This phenomenon was attributed 
to signal amplification. However, since cells may sense 
locally in pseudopods, at least in shallow gradients [ilj ]. 
receptor diffusion may also be responsible for increasing 
the accuracy of sensing. 
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Appendix B: Calculation of £2(0;) 

In this appendix, we calculate Y^ 2 {uj) in Eq. (|29l) . 
Without loss of generality, we set (xo, yo, z o) — (0,0,0) 
and obtain: 



E 2 (w) 



dq±_ f°° qdq iuj/D 2 -q 2 1/D 3 
! Q 2lT J 27T (g2 + ql ) %L q 2 + k_ to 

dq X [°° dq q/D 3 ( k_/D 2 



2n J Q 2n( q *+ q l)-%\ q 2 + tt--%- 
= h(uj) - I 2 (ui), 



where 



, dq±_ f°° qdq l/D 3 k-/D 2 
I l( UJ ) - / TTT 



27T J 2n (qi+ql)- % q 2 + ^ - % 
A qdq i/D 3 k-/D 2 



8tt y/iu/Ds - q 2 q 2 + f- - i£ 



(B2) 



fe_/(87TD 3 ) 



k - ■ I 1 1 \ 



arctan 




A 2 - iuj/D 3 



arctan 



-iu>/D 3 



re- I 1 1 \ / 



and 



. dqj_ f°° qdq 

J 2( W ) = / IT- ' 



o 2tt Jo 2vr -iuj + D 3 {q 2 + q 2 ) 



1 f A qdq 



(B3) 



D 3 Jo 8tt ^/- q 2 + iuj/D 3 

- 1 ( I I^L [~\2 ~^_\ 

~ 8irD 3 \\D~ 3 ~ V ~Da) 

using a cut off to account for the finite size of the receptor. 
Taken together, the integral E 2 (w) is given by: 



Appendix A: Calculation of Ei(oi) 

We devote this Appendix to the calculation of the in- 
tegral £i(<j) in Eq. (|26|) . Ei(u;) is given by: 



Ei(w) 



dq 



1/D 2 



k- 



D-2 



dq q/D2 



2tt q i 



fc- _ ' 

D 2 D 2 

(Al) 

which is divergent. Hence, we introduce a cut-off A ~ 1/s 
to account for the finite dimension s of the receptor. This 
procedure leads to the final result: 



Ei(w) 



dq q/D 2 
2n „2 4. !Ll - iOL 

y T D 2 D 2 



log( 



1 



D2A- 2 



4irD ? 



(A2) 



E 2 (w) = h(uj) - h(w) 
k-/D 2 



A 2 



fc_ 

D 2 



Mih - 777) 



arctan 



arctan 



fc_ 

D 2 



IU1 

D~ 3 



A 2 - 



] 1 



8nD 3 



(B4) 
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